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ABSTRACT 

A  simple  adiabatic  theory  of  particle  modulations  is  derived. 
Assuming  that  the  perturbed  magnetic  field  is  tine  dependent  but 
homogenous  in  space  along  the  dipole  direction,  an  expression  for 
the  change  of  particle  fluxes  in  the  eouatorial  plane  is  derived  in 
terms  of  measurable  quantities.  We  shovr  here  that  the  derived 
expression  agrees  well  with  observed  characteristics  of  trapped 
electron  fluxes  measured  by  the  University  of  Minnesota  electron 
spectrometer  on  the  ATS-1.  In  particular,  quasi-periodic.  modulations 
of  paper  1  are  satisfactorily  explained. 
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l .  Introduct  ion 

In  Part  1,  It  was  shown  that  the  trapped  electron  flux  variations  with 
time  scales  of  a  few  minutes  are  wel 1-organlzed  as  a  single  function  of  1.. 

We  noted  in  Figure  5  of  Part  1,  which  is  reproduced  here  as  Figure  1, 
that  the  enhancement  of  electron  fluxes  is  strongly  dependent  on  the 
equatorial  pitch-angles.  In  particular,  the  largest  changes  are  observed 
near  90°  pitch-angles.  A  conclusion  reached  by  us  was  that  for  the  most 
part  £  50  keV  trapped  electron  fluxes  respond  adiabatieal ly  to  local  magnetic 
field  changes  in  the  evening  magnetosphere. 

When  the  magnetosphere  is  compressed  (inflated)  adiabatieal lv ,  the  particle 
energy  and  the  mirror  points  will  be  altered  (McTlwain,  19bf>;  Kaufmann,  1974; 
Murphy  et.  nl.,  1975).  These  changes  in  particle  energy  and  pitch-angles 
modulate  trapped  particle  fluxes.  If  the  particle  distribution  is  not  uniform 
in  space,  a  detector  fixed  in  space  will  observe  additional  modulations  because 
the  detector  will  be  able  to  scan  a  range  of  1,-shells  (l.ezniak  and  Winckler, 
1970;  Lin  and  Parks,  1974).  The  main  purpose  of  this  paper  is  to  show  that  the 
features  observed  bv  the  electron  spectrometer  reported  in  Part  1  are  consistent 
with  effects  that  arise  when  the  magnetic  field  variations  are  coupled 
adlabatically  to  trapped  electron  fluxes. 

Our  main  interest  is  to  understand  quantitatively  the  variations  of  the 
pitch-angle  distributions  that  occur  within  one  oscillation  period  (see  Figure 
l).  For  this  study,  it  will  be  assumed  that  the  magnetic  oscillations  have 
an  azimuthal  symmetry  in  the  magnetosphere .  This  assumption  is  Justified  if 
the  distance  the  particles  drift  in  one  oscillation  period  is  much  less  than 
the  size  of  the  oscillation  region.  Numer ical ly ,  us ing  the  value  of  the  period 
as  5  minutes  and  100  minutes  for  the  dritt  period,  one  finds  that  the 
oscillation  region  must  be  much  more  extended  than  about  70*  in  azimuth.  This 
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condition  is  satisfied  because  the  data  strongly  suggest  that  the  2-12 
minute  oscillations  are  large-scale. 

If  the  magnetic  variations  have  an  azimuthal  symmetry  at  all  times,  it 
is  known  that  in  the  absence  of  steady  E-field  the  drifting  particles  will 
always  remain  on  the  same  L  shell  (Roederer,  1970).  This  means  that  even 
though  the  2-12  minute  time  scale  is  shorter  than  the  drift  period  of  the 
particles,  the  third  invariant  will  still  be  conserved  during  the  oscillations. 
The  time  scale  of  interest  here  is  much  longer  than  the  cyclotron  or  the 
bounce  period  of  these  particles.  Hence,  the  first  and  the  second  invariants 
are  also  conserved.  In  the  quantitative  formulation  given  below,  it  will  be 
assumed  that  all  of  the  three  invariants  arc  conserved  during  the  oscillations. 

The  approach  taken  here  is  an  extension  of  the  work  previously  published 
by  Mcllwain  (1966)  who  studied  the  adiabatic  behavior  of  trapped  equatoriallv 
mirroring  particles.  Conserving  all  of  the  three  adiabatic  invariants,  our 
formulation  will  include  results  of  particle  energy  and  pitch-angle  changes 
that  occur  for  particles  that  mirror  off  the  equator.  Also,  using  the  l.iou- 
ville  Theorem  and  a  given  initial  spatial  distribution  (which  is  approximated 
from  observations)  we  will  predict  the  ratios  of  particle  fluxes  at  two  different 
times  at  the  same  geographic  location  when  the  magnetic  strength  changes.  The 
theory  is  then  compared  to  observations. 

The  theoretical  formulation  is  presented  in  Section  2.  Application  of 
the  theory  follows  in  Section  3. 

2 .  Theoretical  Forim ilat ion 

2.1  A  Magnetic  Field  Model:  To  describe  the  adiabatic  variations  we  will 
assume  that  the  geomagnetic  dipole  field  is  perturbed  by  b(t).  b(t)  is 
considered  to  be  only  in  the  z-direction,  uniform  in  space,  and  time  varying. 
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The  expression  for  the  total  magnetic  field  taking  into  account  this 
perturbation  is 


B(t)  =  -(2M/r3)cosXr  +  (M/r3)sinXX  +  b(t)z  -  R  r  +  B.X  .  (1) 

r  a 

A  A  A 


M  is  the  dipole  moment  of  the  earth  and  the  coordinates  r,X,  and  7.  are  the 
directional  unit  vectors  defined  in  Figure  2.  The  equation  of  a  magnetic 
field  line  at  any  time  during  the  perturbation  is  obtained  by  integrating 
(rdX/B^)  =  (dr/B^).  Using  F.quation  (1),  one  obtains 

r (t)  -  L  (1  -  br3/2M)cos2X . (2) 

Notice  that  in  the  absence  of  any  perturbation  (b=0) ,  Equation  (2) 
reduces  to  the  familiar  equation  for  the  dipole  field  line. 

Combining  Equations  (1)  and  (2)  one  obtains  B(t,L,X),  the  total  magnetic 


field  strength: 

B(t,L,X)  «(M/L3)(1  +  3sinX)1/2/cos6X 

+  (b/2)  (5  +  3sin2X) /  (1  +  3sin2X)1/2 . .  (3) 

When  X°0,  Equation  (3)  reduces  to 

B(t,L,X»0)-  M/L3  +  5b/ 2 .  (A) 


2.2  Adiabatic  Invariants:  The  expression  of  the  second  adiabatic  invariant 

in  a  time-dependent  magnetic  field  is 

.  Xm  1/2 

K(L,Xm)  -  (M/L)1'2  *(B(t,I,,Xm)  -  B(t,L,X)>  ds . (5) 

where  Xm  is  the  mirror  point  latitude  at  time  t.  This  intergral  must  be  computed 

along  the  perturbed  field  line.  Using  Equation  (2),  one  finds  for  ds 

•  (  ,  2\»i/2  r  1  6cos&X  (1  +  15  sin2X)  1 

ds  -  dXLcosX  (1  +  3sin  X)  1  1 - 2fT+3si"nTX) - -  3  ’  *  ’ 

Here,  6  is  defined  as  bL3/M. 
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Let  Xmo  be  the  particle  mirror  latitude  at  t  -  t  when  there  is  no 
magnetic  perturbation.  Then  the  second  adiabatic  invariant  K(L,Xm)  can  be 
Taylor  expanded  about  the  unperturbed  value  K(L,Xmo): 


K(L,Xm)  =  K(L,Xmo)  +  (Xm  -  Xmo) 


+  6  3K 
Xmo  36 


Xm 

6 


The  conservation  of  the  second  invariant  means  that  K(L,Xm) 
From  Equation  (7),  Xm  and  Xmo  are  related  by 

Xm  *  Xmo  -  6 (3K/36) / (3K/3Xm)  ,  . 


X  mo  ... 
o 

=  K(L,Xmo) 


(7) 


(8) 


It  can  be  shown  that  3K/36  <o.  lienee,  Xm>Xmo  when  6>o.  This  means  that 
the  particle  mirror  latitude  is  lowered  toward  the  earth  when  the  magnetosphere 
is  compressed. 

(Xm-Xmo)/6  is  plotted  in  Figure  3  as  a  function  of  the  equatorial  pitch- 
angle  For  practical  purpose,  the  ratio  (Xm-Xmo)/6  can  be  expressed  with 

an  error  of  less  than  1%  by  the  following 


(Xm-Xmo)/6  =  0.30Sin4Xm  +  0.11Sin2Xm  .  (9) 

Now  we  compute  the  corresponding  change  of  energy  during  the  adiabatic 
compression  using  the  first  adiabatic  invariant.  Let  Eq  be  the  particle 
energy  at  t  and  E  be  the  energy  at  time  t.  Then  it  follows  that 

E/Eq  ■  B(t ,L,Xm) /BtQ,  L,Xmo)  ■  1  +  (l/B)0B/3Xm) (Xm  -  Xmo) 

+  (6/B) (3B/36)  .  (10) 

The  derivatives  are  to  be  evaluated  at  Xm  *«  Xmo  and  6=0.  It  can  be  shown 

using  Equation  (3)  that 

"  OSinXm/cosXm)  *  \  ^>Xm)  . (U) 


t 
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and 


(JL)  (1®) 


cos6Xm 


(5  +  3sin2Am) 
2(1  +  3sin*Xm) 


(12) 


Equation  (10)  can  be  rewritten  as 

E/Eo  -  1  +  6£G(Arao)  +  H(Amo)] . .  (13) 

where  G(Xm)  and  H(Am)  are  defined  as 

Jin 

G(\m)  =  (l/B)(^^)(Xm  -  Xmo)/6 . (14) 

and 

-  III . 05) 


Figure  4  shows  the  plot  of  the  functions  G(Am)  and  H(Am).  Note  that 
K(Am)  is  larger  at  large  pitch-angles  while  G(Am)  Is  larger  at  small  pitch- 
angles.  Phvsicollv,  G(Am)  is  related  to  the  change  of  the  mirror  points  and 
therefore  is  due  to  Fermi  acceleration.  H(Xm)  is  related  to  compressional 
effects  and  therefore  is  due  to  batatron  acceleration.  Note  that  at  30° 
pitch-angle,  the  contribution  from  Fermi  and  betatron  actions  are  equal. 


2.3  Particle  Flux  Modulation:  Let  J(E,o,  t,L)  be  the  differential 

particle  flux  at  energy  F,  and  pitch-angle  a  at  time  t  at  the  magnetic 

coordinate  L.  L  is  related  to  the  position  r  through  Equation  (2).  To 

understand  the  particle  flux  modulation  due  to  adiabatic  effects,  let  y  be 

the  ratio  of  particle  fluxes  at  the  energy  E  and  pitch-angle  a  measured  at  the 

times  t  and  t  at  the  coordinates  L  and  L  : 

o  o 

Y  -  J(F.,a,t,L)/J(E,a,to,l.o) . (16) 
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L  and  Lq  are  the  magnetic  coordinates  at  the  position  r  at  times  t  and 

tQ.  t  cooresponds  to  the  time  when  the  magnetic  perturbation  is  zero. 

Since  our  intention  is  to  understand  observational  features  in  the 

equatorial  plane,  in  the  subsequent  computations  we  will  set  A«0.  From 

Equation  (2),  one  then  notes  that  L  =  r  and  L  ■  r/(l  -  br3/2M). 

o 

Using  Liouville’s  theorem  and  the  fact  that  particles  in  our  model 
remain  on  the  same  L  shell  during  the  perturbation,  we  can  conclude  that 


J(E,a,t,L)  -  <E/Eo)  J(Eo,ao,to,L)  .  .  .  .(17) 

The  relation  between  E  and  E  is  obtained  from  Equation  (13).  The 

o 


relation  between  (ot  )  and  (a  )  can  be  derived  from  the  first  aidabatic 
invariant.  Then, 

sin2a/sinza  »  (E  /E)B(t,L,A  -  0)/B(t  ,L,*  «*  0) . (18) 

o  o  o 

Here  B(t  ,L, A  **  0)  **  M/L3  and  using  Equation  (4),  one  finds  that 


B(t  ,L,  A  =  0)  =  B(tQ,L,A-  0)  +  2.5b. 

Let  us  now  assume  that  J(E,a,to,L)  can  be  represented  as 

J(E,a,to,L)  -  CBa(to,L)E_nsinnia  .  (19) 

Then,  substituting  Equations  (17),  (18),  and  (19)  into  Equation  (16),  one 
obtains 

Y  -  (E  /E)-1~n"m/2  B(t  ,L  )”aB(t  ,L)a_in/2B(t ,L)~m/2 . (20) 

O  O  °  o 


From  Equation  (1),  note  that  B(t,L,A  »  0)  ■  B(tQ,Lo,A  *  0)  +  b(t). 

We  now  use  Equation  (13)  and  the  relationship  between  B(t,L)  and 

B(t  ,L  )  and  B(t  ,L)  to  derive 
o*o  o 

Y  •  1  +  5^(1  +  n  +  m/2)(G(Amo)  +  H(Amo))-  1.5a  -  1.25m} . (21) 


/ 


Page  7 


Note  that  Equation  (21)  reduces  to  Equation  (15)  of  Mcllwain  (1966) 
for  Amo  ■  0  and  m  **  0.  All  quantities  in  Equation  (21)  are  measurable 
experimentally . 

3.  Application  of  the  Adiabatic  Theory 

We  wish  to  compare  Equation  (21)  to  data  obtained  at  synchronous 
altitudes.  Using  Figure  4  of  Paper  1,  wa  obtain  the  value  for  "a"  to  be 
^  4.5.  This  value  was  obtained  for  flux  variations  that  occurred  in  the 
interval  of  time  around  0100  -  0200  U.T.  In  Figure  5,  we  have  plotted 
(AJ/J)/(AB/B)  =■  (6-  1)/S  as  a  function  of  the  equatorial  pitch-angles  of  the 
particles.  The  three  points  shown  in  the  pitch-angle  range  about  70°  -  90° 
are  obtained  from  the  ATS-1  electron  data  and  they  represent  values  of  particles 
whose  energy  is  greater  than  50  keV.  These  points  were  obtained  for  the 
variations  that  occurred  in  the  UT  time  interval  0100  -  0200.  Here  AJ 
represents  the  flux  changes  that  occurred  from  minimum  to  maximum  during 
the  changes  of  the  local  magnetic  field  intensity  AB.  The  values  shown  are 
averages  of  seven  modulation  events  that  occurred  in  this  time  interval, 

(Figure  1  shows  a  typical  example). 

The  experimentally  computed  values  have  been  compared  against  the 
theoretical  curves  obtained  from  Equation  (21).  Since  Figure  1  suggests  that 
m  is  negative,  we  have  used  m  ■  -1  in  our  computation.  With  this  value  of 
m,  we  find  that  the  expeimental  values  agree  best  with  the  theory  for  n 
between  3  and  4  (see  Figure  5). 

The  ATS-1  electron  spectrometer  was  not  designed  to  obtain  accurate 
energy  spectral  estimates  since  the  energy  windows  are  wide.  However,  the 
energy  spectral  index  predicted  by  the  adiabatic  theory  suggests  that  during 


/ 


the  modulation  events  the  index  was  around  3  and  4.  It  is  worthy  of  note 
that  such  index  values  for  trapped  electron  fluxes  are  quite  reasonable  on 
the  basis  of  other  observations  (Mcllwain,  1966;  Lezniak  and  Winckler,  1970; 
Pizzela  and  Frank,  1971). 

The  ATS-1  only  provided  data  in  the  large  pitch-angle  interval.  While 
the  theory  and  the  data  agree  well  here,  we  do  not  have  at  this  time  data 
at  small  pitch-angles  to  verify  whether  the  theory  is  valid  there.  An 
important  point  to  note  is  that  below  about  40°  pitch-angle,  (AJ/J)/ (AB/B) 
becomes  negative  (Figure  5).  This  means  that  the  particle  fluxes  at  small 
pitch-angles  decrease  when  the  magnetosphere  is  compressed,  an  opposite 
behavior  to  the  particles  at  large  pitch-angles.  Physically,  this  effect  can 
be  understood  by  noting  that  the  particle  flux  changes  are  due  to  the  combined 
effects  of  betatron  and  Fermi  accelerations  and  L-shell  modulation. 

4.  Discussion 

In  our  formulation  of  the  adiabatic  modulation  theory,  recall  that  we 
considered  that  the  perturbed  magnetic  field  is  only  in  the  dipole  direction. 
In  reality,  however,  perturbations  are  also  observed  in  the  transverse 
directions  (Barfield  and  McPherron,  1972;  Lin  and  Parks,  1974).  The  fact 
that  the  limited  theory  still  agrees  well  with  the  observations  suggests  that 
the  effects  from  transverse  perturbations  are  small.  Qualitatively,  this 
can  be  understood  by  noting  that  in  the  expression  of  the  first  invariant 
li  *■  E/{(Bq  +  b,, ) 2  +  b^}  the  energy  change  is  proportional  to  b„/Bo  and 

/ B2.  Evidently,  the  second  order  effects  are  negligible. 

*  o 

Additionally,  the  model  presented  here  did  not  include  effects  of 
convection  and  radial  diffusion  in  computing  the  particle  flux  modulations. 


While  this  omission  is  justified  for  studying  particle  modulations  of  a 
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few  minutes  period,  these  effects  must  be  included  if  the  time  scales 
involved  are  comparable  to  convective  and  diffusive  times. 

In  spite  of  the  simplicity  of  our  model,  the  theory  has  provided 
a  means  for  separating  the  adiabatic  and  non-adiabatic  effects  in  the  data, 
making  it  possible  to  study  more  quantitatively  the  contributions  of  non- 
adiabatic  processes.  We  are  currently  in  the  process  of  removing  the 
adiabatic  effects  in  particle  acceleration  events  to  study  in  more  detail 
the  non-adiabatic  contributions. 
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FIGURE  CAPTIONS 


Figure  1:  Behavior  of  pitch-angle  distribution  during  adiabatic 
oscillations.  Large  enhancements  in  fluxes  are  observed  at 
high  pitch-angles. 

Figure  2:  Coordinate  system  used  for  the  geomagnetic  field 
perturbation. 

Figure  3:  Fractional  change  of  the  mirror  latitude  normalized  to 
the  fractional  change  in  the  magnetic  field  strength  as  a 
function  of  the  mirror  latitude  and  equatorial  pitch-angle. 
Corresponding  equatorial  pitch-angle  scale  is  shown  in  the  top 
scale. 

Figure  4:  Fractional  change  of  the  particle  energy  normalized  to 
the  fractional  change  of  the  magnetic  field  intensity  due  to 
betatron  and  Fermi  acceleration.  Note  that  equal  effects  are 
observed  at  about  30°  equatorial  pitch-angle. 

Figure  5:  Fractional  change  of  fluxes  predicted  normalized  to  the 
fractional  change  of  the  magnetic  field  intensity  as  a  function 
of  the  equatorial  pitch-angle.  The  solid  curves  are  derived 
assuming  a  form  of  the  distribution  shown,  "a"  is  a  measure  of 
the  gradient  and  is  derived  from  observations.  Observed  data 
points  from  June  26,  1967  events  are  shown  to  agree  with  the 
theoretical  predications  with  the  parameters  indicated  in  the 
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solid  curves 
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